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HIGH LEVEL EXCURSION SET GEOMETRY FOR 
NON-GAUSSIAN INFINITELY DIVISIBLE RANDOM 

FIELDS 

By Robert J. Adler* Gennady SAMORODNiTSKY*'t 
AND Jonathan E. TAYLOR*'t 

Technion, Cornell and Stanford. 

We consider smooth, infinitely divisible random fields (^X{t), t G 
M), M C R'', with regularly varying Levy measure, and are inter- 
ested in the geometric characteristics of the excursion sets 

Au = {teM : X{i) > u} 

over high levels u. 

For a large class of such random fields we compute the u — > oo 
asymptotic joint distribution of the numbers of critical points, of vari- 
ous types, of X in Au, conditional on Au being non-empty. This allows 
us, for example, to obtain the asymptotic conditional distribution of 
the Euler characteristic of the excursion set. 

In a significant departure from the Gaussian situation, the high 
level excursion sets for these random fields can have quite a compli- 
cated geometry. Whereas in the Gaussian case non-empty excursion 
sets are, with high probability, roughly ellipsoidal, in the more general 
infinitely divisible setting almost any shape is possible. 



1. Introduction. Let {X{t), t G M), where M is a compact set in M*^ 
of a kind to be specified later, be a smooth infinitely divisible random field. 
We shall assume, in a sense that we shall make precise later, that X has 
regularly varying tails. Note that this means that the tails of X are heavier 
than exponential and, in particular, heavier than those of a Gaussian random 
field. Nevertheless, the model we are considering allows both heavy tails (e.g. 
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infinite mean or variance) and light tails, in the sense of the existence of finite 
moments of arbitrary given order. 

We are interested in studying the excursions of the random field over 
levels n > 0, particularly when the level u becomes high. Writing 



for the excursion set of X over the level we shall study the geometric 
characteristics of Au under the condition that it is not empty, i.e. under 
the condition that the level u is, in fact, exceeded. In particular, we shall 
be interested in computing the conditional limit distribution of the Euler 
characteristic of as u — > oo. We refer the reader to [2] for a recent 
detailed exposition of the geometric theory of the excursion sets of smooth 
Gaussian and related random fields, and to [1] for applications of the theory. 

In a significant departure from the well understood Gaussian situation, 
the excursion sets over high levels for the random fields in this paper can have 
quite a complicated geometry. In the Gaussian case excursion sets, unless 
they are empty, tend, with high probability, to contain a single component 
which is almost ellipsoidal in shape, and so have an Euler characteristic 
equal to one. In contrast, the Euler characteristics of the excursion sets in 
our fields can have highly non-degenerate conditional distributions. As a 
consequence, these models are sufficiently fiexible to open the possibility of 
fitting empirically observed excursion sets with widely different geometric 
characteristics. This, more statistical, problem is something we plan to tackle 
in the future. 

The main result of the paper is Theorem 4.1. While it is rather too techni- 
cal to summarise here in full, here is the beginning of a special case. Suppose 
that Nx{i,u) is the number of critical points of X in A^ of index i. Thus, 
if d = 2, Nx{0, u) is the number of local minima of X above the level u in 
the interior of M, Nx{^,u) the number of saddle points and iVx(2,'u) the 
number of local maxima, all above the level u. Then Theorem 4.1 gives an 
explicit expression for the limiting joint distribution 



In fact. Theorem 4.1 goes far beyond this, since it includes not only these 
critical points, but also the critical points of X restricted to the boundary 
of M. The importance of this result lies in the fact that Morse theory shows 
how to use the full collection of these critical points to describe much of the 
geometry of A^, whether this geometry be algebraic, integral, or differential. 



(1) 



Au = Au{X, M) = {teM : X{t) > u} 
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Furthermore, Theorem 4.1 can also be exploited to describe a very simple 
stochastic model for high level excursion sets, as well as to develop a simple 
algorithm for simulating them. 

The remainder of the paper begins in Section 2 with a description of the 
parameter spaces M with which we work, and is then organised as follows: 
In Section 3 we define our model, discuss the smoothness assumptions we 
are imposing, as well as those related to the regular variation of the tails. 
Section 4 contains the main result of the paper, on the joint distribution of 
the numbers of high level critical points of infinitely divisible random fields. 
This is followed with one of its main consequences, the distribution of the 
Euler characteristic of high level excursion sets, in Section 5. In Section 6 
we introduce a class of moving average infinitely divisible random fields and 
derive conditions under which the main result of the Section 4 applies to 
them. We also provide examples to show that, by choosing appropriately 
the parameters of the model, one can make the geometric structure of the 
high level excursion sets either 'Gaussian-like' or 'non-Gaussian-like'. Fi- 
nally, Section 7 contains the proof of the main theorem for the special case 
of M the unit cube in W^, while Section 8 discusses how to extend this to 
more general parameter spaces. 

Throughout the paper C stands for finite positive constants whose precise 
value is irrelevant and which may change from line to line. 

2. The parameter space M. Throughout this paper we shall assume 
that M is a compact, C^, Whitney stratified manifold embedded in M"^. 
These are basically sets that can be partitioned into the disjoint union of 
manifolds, so that we can write M as the disjoint union 



where each stratum, djM, < j < dim(M), is itself a disjoint union of 
a number of bounded j-dimensional manifolds. We shall also assume that 
each of these manifolds has bounded curvature. The "Whitney" part of the 
definition relates to rules as to how the djM are glued together. 

The simplest example of such a stratified manifold would be a simple 
manifold, in which there would be only a single term, of maximal dimension 
dim(M), in (3). A domain M C M*^ is not much more complicated, 
and could be stratified into its interior, of dimension d, and boundary, of 
dimension d — 1. 

A familiar example of a stratified manifold which has components of all 
dimensions is the d-dimensional cube 1^ = [—1, l]"^- In this case d^M would 



dim M 



(3) 




i=o 
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be the interior of the cube, dd~iM the union of the interiors of the 2d 'sides' 
of dimension (d — 1), and so on, with BqM being made up of the 2'^ vertices. 

For further details on Whitney stratified manifolds see, for example, Chap- 
ter 8 of [2]. However, if you want to avoid this level of generality, you will 
lose little (beyond many applications) by taking M = throughout the 
paper. In fact, the main proofs in Section 7 treat only this case in full, with 
only a brief description, in Section 8 of how to extend the result to general 
M. 

3. Smooth infinitely divisible random fields and regular varia- 
tion. In this section we shally define the random fields of interest to us, 
describe their distributional structure, and then specify the smoothness as- 
sumptions necessary for studying the geometry of their excursion sets. 

A reader familiar with the theory of infinitely divisible processes will note 
that the route we take goes back to first principles to some extent. (For 
example, it would be more standard, nowadays, to start with the function 
space Levy measure Ax of Section 3.3 rather than invest a couple of pages in 
defining it.) The need for this, as should become clear below, is to be able to 
carefully define random fields, along with their first and second order partial 
derivatives, on a common probability space. 

3.1. Probabilistic structure of infinitely divisible random fields. As a first 
step, we shall need to define our random fields on a region slightly larger 
than the basic parameter space M, and so we take M be a bounded open 
set in R'^, with M cM. 

We now consider infinitely divisible random fields of the form 

(4) X{t) = f f{s;t)fi{ds), teM, 

Js 

where {S,S) is a measurable space and /i is an infinitely divisible random 
measure on S with characteristics defined below. (We refer you to [11] for 
more information on infinitely divisible random measures and stochastic in- 
tegrals with respect to these measures. See also the forthcoming monograph 
[17] for a more complete account.) 

The infinitely divisible random measure fi, which we shall define in a 
moment, is characterised by its 'generating triple' (7, -F, /3). Here 7 is a u- 
finite measure on (S,S), and plays the role of the variance measure for the 
Gaussian part of fi. More important for us is the Levy measure F, which 
is a (j-finite measure on S" x (M \ {0}), equipped with the product a-field. 
Finally, /3 is a signed measure on (5,5), which plays the role of the shift 
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measure for //. Denote by So the collection of sets B in S for which 

^{B) + m\{B)+ f [[xf F(i?,dx)<(X), 
Jr\{o} 

where 




X if |x| < 1, 

sign(2;) otherwise, 



and is the total variation norm of /?. 

With all elements of the triple defined, we can now define the infinitely 
divisible random measure {^{B), B £ Sq) as a stochastic process for which, 
for every sequence of disjoint 5o-sets Bi,B2, ■ ■ ., the random variables /^(-Bi), 
IJ.{B2), . . . are independent (that is, fj, is independently scattered) and if, in 
addition, [j^Bn € Sq, then n{[J^Bn) = J2nf^iBn) a.s. (that is, fi is a- 
additive). Finally, for every B G Sq, /i(-B) is an infinitely divisible random 
variable with characteristic function given by 

= exp I - ^7(5)02 + [ {e'^'' -l-ieix]])F{B,dx)+i9P{B)}, 
^ ^ Jr\{o} > 

for ^ € M. The monograph [19] can be consulted for information on infinitely 
divisible random variables. 

We shall assume (without loss of generality) that the Levy measure F has 
the form 

(5) F{A) = [ p{s;As)m{ds), 

Js 

for each measurable Ac Sx (M\{0}) , where Ag = {x e M\{0} : (s, x) G A} 
is the s-section of the set A. In (5), m is a u-finite measure on {S,S) (the 
control measure of fi), and the measures {p{s; •)) (the local Levy measures) 
form a family of Levy measures on M such that for every Borel set C C 
M\{0}, s — > p(s; C) is a measurable function on S. We can, and shall, choose 
the control measure m in (5) in such a way that \\(3\\ is absolutely continuous 
with respect to m, and define the Radon-Nikodym derivative b = dp/dm. 
The local Levy measures p, which, intuitively, control the Poisson structure 
of the random measure fi around different points of the space S, will play a 
central role in all that follows. 

Note that while it is possible, and common, to choose m in with the 
added feature that 7 is also absolutely continuous with respect to m, and 
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that p(s;M \ {0}) > on a set of s € 5 of full measure m, we shall not 
require this and so shall not do so. 

Finally, we assume that the kernel s € /S*, t G M in (4) is deter- 

ministic and real, such that, for every t G M, the mapping f{-;t) : S ^ M 
is measurable, and that the following three inequalities hold: 



(6) / fis;tyj{ds) < oo, 

Js 

(7) /_ f Jxf{s;t)fF{ds,dx) < oo, 
and 

(8) / bis)fis; t) + / {[[xf{s; i)]] - [[x]]f{s; t)) p{s; dx) m{ds) < oo , 



These conditions guarantee that the random field [X[t), t G M) in (4) is 
well defined. 

A particularly simple, but rather useful, example of this setup is studied in 
Section 6 below, when X is a moving average random field. In this example 
both 7 and /3 components of the generating triple vanish, so, in particular, 
the random field has no Gaussian component. Furthermore, S = M^, the 
control measure m is Lebesgue, and the local Levy measures p{s, •) are inde- 
pendent of s. Finally, the kernel function / is of the form /(s, t) = g(s + t) 
for some suitable g, and so the random field is given by 

(9) X{t)= I g{s + t)fi{ds), tGMcM'^. 

The random measure p has, in this case, the stationarity property n{A) = 
li{t + A) for all Borel ^ of a finite Lebesgue measure and t G M'^, which 
immediately implies that a moving average random field is stationary. An 
impatient reader, who already wants to see results without wading through 
technicalities, might want to now skip directly to Section 6.2 to see what 
our results have to say for moving averages. 

Returning to the model (4), note that it has been defined in considerable 
generality, so as to allow for as wide a range of applications as possible. For 
example, we retain the Gaussian component of the random field X. However, 
the tail assumptions imposed below will have the effect of ensuring that the 
Gaussian component will not play a role in the geometric structure of high 
level excursion sets. 
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3.2. Regularity properties. We shall require that the sample paths of X 
satisfy a number of regularity properties for the theory we are developing 
to hold. The main assumption will be that the paths of X are a.s. C^, for 
which good sufficient conditions exist. The secondary assumptions require 
a little more regularity, centered around the notion of Morse functions. For 
more details see Chapter 9 in [2]. 

Recall that a function / : M ^ R is called a Morse function on the 
stratified manifold M if it satisfies the following two conditions on each 
stratum d^M, k = 0, . . . , dim(M). 

(i) f\d^.M it is non-degenerate on dkM, in the sense that the determinant 
of the Hessian of /la^A/ ^-t its critical points does not vanish. 

(ii) The restriction of / to dkM = \Jj=odjM has no critical points on 

Here is our first, and henceforth ubiquitous, assumption. 

Assumption 3.1. On an event of probability 1, the random field X has 
C"^ sample paths on M and is a Morse function on M . 

Sufficient conditions for Assumption 3.1 to hold are not hard to come by. 
As far as the assumption is concerned, it suffices to treat the Gaussian 
and non-Gaussian components of X separately. For the Gaussian part, there 
is a rich theory of necessary and sufficient conditions for continuity, which are 
easy, at least as far as sufficiency is concerned, to translate into conditions 
that ensure the C"^ property also holds. For details, see Section 1.4.2 of [2]. 

Necessary and sufficient conditions for the assumption on the non- 
Gaussian component are not known, but a number of sufficient conditions 
exist. It is not our goal in this paper to develop the best possible conditions 
of this sort, so we restrict ourselves to one situation that covers, nonetheless, 
a wide range of random fields. Specifically, we shall assume that the 7 and 
13 components in the generating triple of the infinitely divisible random 
measure M vanish, and that the local Levy measures p in (5) are symmetric; 
i.e. p{s] —A) = A) for each s € S and each Borel A S M \ {0}. That is, 
M is a symmetric infinitely divisible random measure without a Gaussian 
component. 

The following result gives sufficient conditions for a symmetric infinitely 
divisible random field without a Gaussian component to have sample func- 
tions in C^. (The conditions are also necessary after a slight tightening of 
the assumptions on the null sets involved, cf. Theorem 5.1 of [4].) 
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Theorem 3.2. For a symmetric random field of the form (4), with M an 
infinitely divisible random measure without a Gaussian component, suppose 
that the kernel f : S x M ^ M is (product) -measurable. Assume that for 
every s € S outside of set of zero m-measure the function f{s; •) : M — >■ M 
is C^. Furthermore, assume that the partial derivatives 

df d'^f 
fi{s;t) = —{s;t), i = l,...,d, fij{s]t) = -^^-^{s;t), i,j = l,...,d 

satisfy the following conditions. 

(i) The integrability condition (7) holds when the kernel f{s;t) there is 

replaced by any of the fi{s;t) or fij{s]t). 
(a) The random fields 



(10) X,j{t) = fij{s;t) fiids), t€M, 

s 

i,j = l,...,d, are all sample continuous. 

Then the random field {X{t), t G M) has (a version with) sample functions 
in C72. 

Proof. For i = 1, . . . ,d, define the random fields 

Mt) = [ fi{s]t)n{ds), t G M, 
Js 

where t = {ti, . . . ,td) G M. Furthermore, defining for j = 1, . . . , d 
ej{t) =inf |s : {ti,...,tj-i,s,tj+i,...,td) G m| , 

note that 

Xi{t) = Xi(ti, . . . ,tj-l,ij{t),tj^l, . . . ,tci) 

+ / Xij(ti, . . . ,tj-i, s,tj+i, . . . ,td) ds. 

(cf. Theorem 5.1 of [4].) Since each Xij was assumed to be sample continuous, 
it follows that the Xj are C^. 

Repeating the argument gives that 

X{t) = X{ti, . . . ,ti-i,ij{t),ti+i, . . . ,td) 

+ / Xi(ti, . . . ,tj-i,s,tj+i, . . . ,td) ds , 

which shows that X is C"^. □ 
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Thus, in searching for sufficient conditions for the a.s. second order differ- 
entiabihty of X, it suffices to estabhsh the continuity of the random fields 
of (10). While there are no known necessary and sufficient conditions for 
sample continuity of general infinitely divisible random fields, various suffi- 
cient conditions are available. See, for example. Chapter 10 of [18] for the 
special case of stable random fields, or [9] for some other classes of infinitely 
divisible random fields. 

This is as far as we shall go at the moment discussing the issue of dif- 
ferentiability in Assumption 3.1. Conditions sufficient for X to be a Morse 
function, also required in this assumption, are, in principle, available as well. 
For example, it follows from the arguments of Section 11.3 of [2] (cf. Theo- 
rem 11.3.1 there) that a field X will also be, Morse function on 
the unit cube if the following two conditions are satisfied, for each face J 
of Id, and for all t G J. 

(i) The marginal densities pt{x) of VX| j(t) are continuous at 0, uniformly 
in t. 

(ii) The conditional densities pt{z\x) of Z = det V^X| j(t) given VX| j(t) = 
X are continuous in {z, x) in a neighbourhood of 0, uniformly in t. 

Related conditions are given in Theorem 11.3.4 there for the case of M a 
stratified manifold. 

It does not seem to be trivial to translate the above conditions into general 
conditions on the kernel / and the triple (7, F, /3), and we shall not attempt 
to do so in this paper. On the other hand, given a specific kernel and triple, 
they are generally not too hard to check. In the purely Gaussian case, simple 
sufficient conditions are provided by Corollary 11.3.2 of [2], but it is the more 
involved infinitely divisible case that is at the heart of the current paper. 
If the latter random field is, actually a so-called type-G random field (see 
[15]) (symmetric a-stable random fields, < a < 2 are a special case of 
type-G random fields), then these fields can be represented as mixtures of 
centered Gaussian random fields, and Corollary 11.3.2 in [2] may be helpful 
once again. 

We close this section with a remark and a further assumption. 

Remark 3.3. Unless X is Gaussian, Assumption 3.1 implies that it is 
possible to modify the kernel / in (4), without changing the finite dimen- 
sional distributions of X, in such a way that f{s, •) is G^ for every s £ S; 
see Theorem 4 of [13]. For simplicity we shall therefore assume throughout 
that / has such G^ sections. This ensures, in particular, measurability of 
functions of the type sup^g^ |/(s, i)|, s € 5, which we shall take as given in 
what follows. 



10 



ABLER, SAMORODNITSKY AND TAYLOR 



Assumption 3.4. The kernel f{s,t), s e S, t e M, along with its 
first and second order spatial partial derivatives fi and fij are (uniformly) 
bounded and, for for every s G 5, the function f{s, •) is a Morse function 
on M. 

3.3. The function space Levy measure. Although the infinitely divisible 
random fields we are studying in this paper were constructed above via 
stochastic integrals (4) and, as such, are characterised by the triple (7, F, (5) 
of the random measure ji and the kernel /, in what follows the most im- 
portant characteristic of the infinitely divisible random field (4) will be its 

function space Levy measure. This is a measure on the cylinder sets of M^, 
related to the parameters in the integral representation of the field by the 
formula 

(11) \x = Fo TJ^ , 

where F is the Levy measure of the infinitely divisible random measure /i 
and Tf : 5 x (R \ {0}) ^ M^^ is given by 

(12) Tf{s,x)=xf{s,-), xeM\{0}, 

cf. [11]. Thus the finite dimensional distributions of X are given via the joint 
characteristic function 



(13) E|exp{i5]^7jX(t,)}| =exp|-Q(7i,. 



.7fej 



3 

k 



7=1 7=1 



\xidx) + iL(7i, . . . ,7fc) 



for A; > 1, ti, . . . , ffc G M, and real numbers 71, ... , 7^, where Q is a quadratic 
function (corresponding to the Gaussian part of X), and L is a linear func- 
tion (corresponding to the shift). Their exact forms are not important for 
us at the moment. 

Note that the Levy measures of the first and second order partial deriva- 
tives Xi and Xij are similarly (cf. Theorem 5.1, [4]) given by 

(14) X^^=FoT7\ Ax„=FoT7i, i,j = l,...,d. 



3.4. Regular variation. We now turn to the final set of technical assump- 
tions on our infinitely divisible random fields, these being related to the reg- 
ular variation of their Levy measures, and which we formulate in terms of 
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the local Levy measures of (5). These are our final set of assumptions, and 
our main results hinge on them. 

Recall that a function / is regularly varying at infinity, with exponent a, 

if 

(15) lim = A° for all A > 0. 

x->oo j[x) 

Assumption 3.5. There exists a H : {0,oo) ^ (0, oo) that is regularly 
varying at infinity with exponent —a, a > 0, and non-negative measurable 
functions w+ and on S such that 

,. P{s;{u,oo)) p(s;(-oo,-n)) 

(16) hm — = w+[s), Inn -— = W-[s) 

n^oo H[u) u^oo H{u) 

for all s S. Furthermore, the convergence is uniform in the sense there is 
uq > such that, for all u > uq and all s G S, 

p(s;(n,oo)) p{s;{-oo,-u)) 

-— <2w+is), — <2w-{s). 

H[u) 

The following simple lemma relates Assumption 3.5 to the corresponding 
behaviour of the Levy measure Ax on a set of crucial importance to us. 
We adopt the standard notation a+ = max(a,0) and a_ = (— a)+ for the 
positive and negative parts of a real. 

Lemma 3.6. Let Assumption 3.5 hold. 

(i) Assume that the kernel f{s,t), t £ M is uniformly (in s & S) bounded, 
and that for some e > 0, 

(17) / (u;+(s) + u;_(s))sup|/(s,t)|"-'m(ds) <oo. 

JS tGM 

Then 

(18) lim ^{^'^PteMX{t)>u} 



H{u) 

\x\^g : supigjij c/(t) > n| 



lim 



Wj^{s) sup /(s, t)" + w^{s) sup /(s, t)t m{ds), 
teM teM 
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where M can be replaced with M throughout. Furthermore, 

{sup^^~|X(t)| >n} 



(19) 



u— >oo 



lim 

n— >oo 



H{u) 

Xx[g ■■ sup^gj,7 \g{'t)\ > uj 



{w+{s) + w-{s)) sup |/(s, m{ds). 

(ii) Assume that the first order partial derivatives fi{s,t), t G M, i 
1, . . . ,d are uniformly (in s & S) bounded, and that for some e > 0, 



(20) 
Then 

(21) 



{w+{s) + w-{s)) sup \ fi{s, t)\" ^ m{ds) < oo . 



lim 

ti— »oo 



Ax, {9 : sup^gjv7 \9{t)\ > 



H{u) 

{w+{s) + w_{s)) sup \ fi{s,t)\°m{ds). 

t€M 



(iii) Assume that the second order partial derivatives fij{s,t), t € M, 
i,j = 1, . . . ,d are uniformly (in s ^ S) bounded, and that for some e > 0, 



(22) 
Then 

(23) 



[w-^{s) + w-{s)) sup \fij{s, t)\^ ^ m{ds) < 00 . 



lim 



Ax,, {5: sup^6mI5 



(t)| > u} 



H{u) 

{w+{s) +w_{s)) sup \fij{s,t)\°' m{ds). 

t£M 



Proof. The first equality in (18) follows from the second equality there 
by Theorem 2.1 in [16]. As for the second equality in (18), it follows from 
(11) and (5) that 



Axis' : sup 5(t) > u| = / P\s; 



— U 



00 



P s; -00 



sup^^^ f{s,t). 



m{ds) 
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Using the uniform boundedness of the kernel and Potter's bounds (cf. [10] 
or [3], Theorem 1.5.6) we see that for any e > there is C > such that for 
aU n > 1, 



The Hmit (18) now follows from Assumption 3.5 via (17), regular variation, 
and dominated convergence. The proof of (19) is identical, as are the proofs 



Remark 3.7. The assumption of uniform boundedness of the kernel / in 
(4) and its partial derivatives will be kept throughout the paper (it is already 
a part of Assumption 3.4), but the only place it is used is in Lemma 3.6. It 
is not difficult to see that this assumption can be removed at the expense of 
appropriate assumptions on the behaviour near the origin of the local Levy 
measures in (5) and of slightly modifying the integrability condition (17). 
Given that this paper is already rather heavy on notation, we shall continue 
to work with uniform integrability, which helps keep things comparatively 
tidy. Note that it is also clear that, for the purpose of proving (18) alone, 
the integrability assumption (17) could be relaxed. 

4. Limiting distributions for critical points. Our initial aim, as 
described in the Introduction, was to obtain information about the distri- 
bution of the Euler characteristic of the excursion sets of (1). As is known 
from Morse critical point theory, Euler characteristics of excursion sets are 
closely related to the critical points above fixed levels. We shall describe this 
connection in the following section and, for the moment, concentrate on the 
critical points of X, which are also of intrinsic interest. 

Recall the partition (3) of the stratified manifold M into collections dkM 
of manifolds of dimension k. Let J denote one such manifold, of dimension 
< A; < d = dim(M). _ 

Now let 5 be a function on M, and and for z = 0, 1, . . . ,dim(J), let 
Cg{J;i) be the set of critical points of index i of g^j. These are the points 
for which 'Vg{t) is normal to the tangent plane to J at t, and for which the 




and 




of (ii) and (iii). 



□ 



14 



ABLER, SAMORODNITSKY AND TAYLOR 



index of the Hessian of 5| j, computed with respect to any orthonormal basis 
for the tangent plane, when considered as a matrix, has index i. (Recah 
that the index of a matrix is the number of its negative eigenvalues.) For a 
full discussion of the critical points of functions on stratified manifolds see 
Chapter 9 of [2]. Let 

Ng{J;i) = Card{Cg{J;i)), 

and, for real u, 

Ng{J\i : u) = Cavd{Cg{J;i) (1 {t : g(t) > u}), 

be the the overall number of the critical points of different types of g, and 
the number of these critical points above the level u, correspondingly. Since 
g is a Morse function, it is standard fare that all of the above numbers are 
finite, (e.g. [2]) 

Just a little more notation is required for the main theorem. Let / be 
the kernel in the integral representation (4) of an infinitely divisible random 
field. For A; = 0, 1, . . . , d, a manifold J and i = 0, 1, . . . , dim( J), let 

Ci{J;s) = iV/(^..)(J;i) 

be the number of the critical points of the s-section of / of the appropriate 
type, well defined since by Assumption 3.4 the sections are Morse functions. 

Furthermore, let {ti{J; i; s), I = 1, . . . , Ci{J; s)) be an enumeration of these 
critical points, and, for 1 < m < Ci{J;s) let 

f[m] ^^(•5)) f[m] ^(^)' 

be, correspondingly, the m-th largest of the positive parts (/(s; ti{J; i; s))) ^, 
1 = 1,..., Ci{J; s), and the m-th largest of the negative parts (/(s; t/( J; i; s))) _ 

/ = 1, . . . , Cj( J; s). (Both of these quantities are set to be equal zero if 
m > Ci{J;s)). 

Finally, extend these definitions to m = by setting 

/(f +)(s) = sup(/(.;t))^, /if -^(s) = sup(/(.;t))_. 

The following theorem is the main result of this paper. It describes the 
limiting, conditional, joint distribution of the number of critical points of 
all possible types of a infinitely divisible random field over the level u, as 
ti ^ oo, given that the random field actually exceeds level u at some point. 
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Theorem 4.1. Let {X{t), t € M) be an infinitely divisible random field 
with representation (4), satisfying Assumptions 3.1, 3.4 and 3.5. Assume, 
furthermore, that (17) holds for some e > 0. Then, for any collection J of 
manifolds in the various strata d^M , k £ {0, 1, . . . , d}, of M , and any collec- 
tion of non-negative integers {n(J; i) = 0,1, . . . , i = 0,1, . . . , dim(J), J € 

J}, 



(24) 



F{Nx{J;i : u) > n{J;i), J ej,i = 0, 1, . . . , dim( J) supX(t) > u 

teM 



Is 


w+{s 


) (mill j,i (s)) + (s) (mill j,i /f^^('_J. . J ^ 


is)) 


m{ds) 


Is 


w+{s) suptgjvf fis, t)% + w- (s) sup^gAf f{s, t)- 


m{ds) 



as u 



oo. 



The proof of the theorem appears in Section 7. 

Remark 4.2. While the structure of (24) might be rather forbidding at 
first sight, its meaning is actually rather simple. The main point of Theorem 
4.1 is that, once the random field reaches a high level, its behaviour above 
that level is very similar to that of the much simpler random field. 



(25) 



Z{t) = Vf{W,t), tGM, 



where {V, W) £ (M \ {0}) x 5 is a random pair, the joint law of which is the 
finite restriction of the Levy measure F to the set 

{{x,.s) G (M\ {0}) X S : sup|a;/(s;t)| > l} , 

teM 

normalized to be a probability measure on that set. 

Remark 4.3. In fact, one can go much further than in the previous 
remark, and interpret the limit (24) as showing that limiting conditional 
joint distribution of critical points is a mixture distribution, that can be 
described as follows. Set 



H 



w+{s) sup f{s, t)" + W-{s) sup f{s, t)" m{ds) . 
teM teM 



1. Select a random point W £ S according to the probability law rj on S 
with 



dm 



H- 



w+{s) sup f{s,t)% + w-{s) sup/(s,t)" 
teM teM 



seS. 
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2. Given W = s, select a random value / € {—1,1} according to the law 

P(I = l\W = s) = w+{s)supt^Mfis,t)% 

w+{s) supigA/ f{s, t)% + w^{s) supigjv^ f{s, t)'^ ■ 

3. Let be a random variable independent of W and /, with P{V^ < 
x} = for < X < 1. Then the numbers of critical points {Nx{J',i '■ 
u), J ^ J ■, given that supX > u, have, as n — > oo, the same distribu- 
tion as the numbers of critical points of the random field 



above the level VJ^, if / = 1, and the numbers of critical points of the 
random field 

above the level Vn if / = —1. 



Remark 4.4. While Theorem 4.1 counts critical points classified by 
their indices, there are also other properties of critical points that are of 
topological importance. For example, in [2] considerable emphasis was laid 
on the so-called 'extended outward critical points', these being the critical 
points t £ M for which V f{t) G Nt{M), where Nt{M) is the normal cone 
of M at t. 

Extended outward critical points play a major role in Morse theory, in 
terms of defining the Euler characteristics of excursion sets. It will be easy 
to see from the proof of Theorem 4.1 that its the statement remains true if 
one replaces critical points by extended outward critical points. This will be 
used in certain applications of Theorem 4.1 below. 

5. The Euler characteristic of excursion sets. One application of 
Theorem 4.1 is to the Euler characteristic ip{Au) of the Gxcursion set 
over a high level u. We shall not define the Euler characteristic here, but 
rather send you to [2] for details. The Euler characteristic of an excursion 
set of a Morse function is equal to the alternating sum of the numbers of 
extended outward critical points of the function over the level. This leads to 
the following result, an immediate corollary of Theorem 4.1, (9.4.1) in [2], 
and Remarks 4.3 and 4.4 above. 
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Corollary 5.1. Under the conditions of Theorem 4-1, the conditional 
distribution of the Euler characteristic of the excursion set of an infinitely di- 
visible random field computed with its limiting conditional distribution given 
that the level is exceeded, is given by the mixture of the Euler character- 
istics of the random fields (26) and (27), with the mixing distribution as 
described in Remark 4-3. In particular, the expected Euler characteristic of 
the excursion set of the limiting (conditional) random field is given by 



w+{s) sup f{s, t)1E{C+{s)} + w-{s) sup /(s, t)^E{C_(s)} m{ds) 



Here, for s G 5, C±{s) is the Euler characteristic of the excursion set of the 
field {f{s,t)±/ sup^gj,/ /(s, r)±, t G M) above the level Va- 

6. An example: Moving average fields. The power and variety of 
the results of the previous two sections can aheady be seen in a relatively 
simple but application rich class of random fields, the moving average fields 
that were already introduced at (9). Recall that these are stationary fields 
with representation 

(28) X{t)=( g{s + t)fi{ds), t G M'^. 

In this section we shall consider a moving average random field on M = 
Id = [—1, l]'^, and the ambient manifold is simply M = W^. 

Our basic assumptions, that will hold throughout this section, are 

(i) The function g is on M"^ and satisfies (7) and (8). 

(ii) fj, is an infinitely divisible random measure on M*^, for which the Gaus- 
sian and shift components in the generating triple, 7 and /3, vanish. 

(iii) The control measure m in (5) is d-dimensional Lebesgue measure. 

(iv) The local Levy measures p{s, •) = p{-) are independent of s G M'^. 

By choosing different kernels g we shall see that quite different types of 
high level excursion sets arise, as opposed to the Gaussian case, in which 
ellipsoidal sets are, with high probability, ubiquitous, 

6.1. Checking the condtions of Theorem 4-1 for type G moving averages. 
In this subsection we exhibit a broad family of moving average random fields 
(28) for which we shall verify the conditions required by the main result of 
Section 4. These are the so-called type G random fields. We emphasize that 
the applicability of our main results is not restricted to type G random 
fields. For the latter we can use standard tools to check the assumptions of 
Theorem 4.1, which is why they are presented here. The main result of this 
subsection is 
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Theorem 6.1. A moving average infinitely divisible random field X sat- 
isfying Conditions 6.2 and 6.3 below also satisfies the assumptions of Theo- 
rem 4-1- 

Condition 6.2. The local Levy measure p is a symmetric measure of 
the form 

(29) piB)=E{po{Z~'B)}, 

where B is a Borel set, Z is a standard normal random variable, and po is a 
symmetric Levy measure on M. Furthermore, the function po{{u, oo)), u > 
is regularly varying at infinity with exponent —a, a > 1, and there is (3 ^ 
[1, 2) such that 

po{{u, oo)) < au~'^, 
for all < u < 1, for some < a < oo. 

In fact, for any Levy measure po on M, (29) defines a Levy measure on 
M; see e.g. Proposition 2.2 in [8]. Furthermore, it is simple to check that 
the behaviour of the measures p and po are similar at zero and infinity. 
Specifically, 

(30) 

«-»°opo(('U,Oo)) 

and 

(31) p{{u, oo)) < jmaxdZl'^, l)aii"'^} 

for < M < 1. In particular, a moving average infinitely divisible random 
field satisfying Condition 6.2 automatically also satisfies Assumption 3.5. It 
suffices to choose H{u) = po{{u, oo)), u > and w+{s) = w^{s) = E {Z^}. 

It is condition 6.2 that makes the random field a "type G random field". 
It implies that the random field X can be represented as a ceratin mixture 
of stationary Gaussian fields, cf. [8]. Under the conditions we impose, each 
one of the latter is a.s. a Morse function, which will tell us that the moving 
average itself has sample functions which are, with probability 1, Morse 
functions. 

If it is known from other considerations that the sample functions of a 
specific infinitely divisible random field are, with probability 1, Morse func- 
tions, then (29) is not needed, and only the assumptions on the behaviour of 
the tails of the Levy measure p{{u, oo)) as u — > or u — > oo are required. In 
the present form of Condition 6.2 these assumptions become the conclusions 

(30) and (31) from the corresponding assumptions on the Levy measure pQ. 
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Condition 6.3. The kernel g is in , and its restriction to any bounded 
hypercube is a Morse function. Assume that the first and the second deriva- 
tives gi, i = 1, . . . ,d and gij, i,j = 1, . . . ,d satisfy (7) and (8). Assume, fur- 
ther, that for almost every s there is no subspace of dimension stricly 
less than {d? + 3d)/2 to which the vectors {gi{s), i = 1, . . . ,d, gij{s), i,j = 
1, . . . ,d, i < j) belong. 

Finally, assume that the function 

Tg{s) = sup \g{s + t)\, seM.'^ 
te[-i,i]'* 

satisfies Tg G L°'~^{W^), while the function 

fg{,s)= max \gij{,s)\-\- sup \gijkis + 1)\, s eR'^ 

te[-i,i]'*,*,i,fcei,...,d 

satisfies fg G L°~^(IR'') n L^{R'^) for some e > and for the a and (3 for 
which Condition 6.2 holds. 

Since these are assumptions on the kernel g in the integral representation 
(28) of the random field, and the kernel is often explicitly given, the above 
conditions are, generally, easy to apply. See the examples below. 

Clearly, a moving average infinitely divisible random field satisfying Con- 
dition 6.3 will also satisfy Assumption 3.4. It also satisfies (17). Theorem 
6.1 is then an immediate consequence of the following two lemmas. 

Lemma 6.4. A moving average infinitely divisible random field satisfying 
Condition 6.3 has sample paths in C^. 

Proof. The only assumption of Theorem 3.2 that is not stated explicitly 
above is the sample continuity of the second order derivative random fields. 
To establish that we use Theorem 3.1 in [9]. The conditions of the latter 
theorem are easily seen to follow from the simple facts that, firstly, for any 
function h on W^, 

hit) - hir) , , 
sup --f \r^<d sup \hi{t)\, 

\\t-r\\ t€ld,i=l,-,d 



and, secondly, that the metric entropy condition in Remark 2.1 in [9] trivially 
holds for any q > 2 for the hypercube Id = [—1, 1]'^ and Euclidian distance 
on that hypercube. □ 
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To complete the proof of Theorem 6.1 we need to check that a moving av- 
erage satisfying Conditions 6.2 and 6.3 has sample functions that are, with 
probability 1, Morse functions. As mentioned above, we shall accomplish 
this by representing the random field X as a mixture of zero mean Gaus- 
sian random fields, each one of which will have, with probability 1, sample 
functions that are Morse functions. 

Lemma 6.5. A moving average infinitely divisible random field satisfying 
Conditions 6.2 and 6.3 has sample functions that are, with probability 1, 
Morse functions. 

Proof. Let u and v be probability measures on M and absolutely con- 
tinuous with respect to the Levy measure pQ and to d-dimensional Lebesgue 
measure A^, respectively. Let 

dv , ^ ^ , , du , . 

^P{x) = —{x), x£R, ip{s) = —{s), seM'^. 
dpo dXd 

Then the random field X has a representation as an infinite sum of the form 

oo 

(32) X{t) = Y^ZkVkg{t + Hk)l{mMHk)Tu<l), t € M^ 

k=l 

where (Zi, Z2, . . .) are i.i.d. standard normal random variables, (Vi, V2, . . .) 
are i.i.d. random variables with a common law v, (Hi, H2, . . .) are i.i.d. 
random vectors in with a common law u, and (Fi, . . .) are the points 
of a unit rate Poisson process on (0, 00). All four sequences are independent. 
See Section 5 in [8] for details. Furthermore, by Lemma 6.4 and Theorem 
3.2, the first and second order partial derivatives of X are also moving 
average random fields, and so have a corresponding series representation. In 
particular, 

00 

(33) X^it) = J2 ZkVk9i{t + Hk)l{^{VkMHk)Tk < 1), t G M'^ 

k=i 

for each i = 1, . . . ,d and 

00 

(34) X,,{t) = ZkVkg^J{t + i/fe)l(V'(Ffc)(^(i/fc)Ffe < 1), f G M'^ 

k=l 

for each i,j = l,...,d. 

We may assume, without loss of generality, that the standard Gaussian 
sequence (Zi,Z2,...) is defined on a probability space (f^i, JFi,Fi), and 



EXCURSION SET GEOMETRY 



21 



the remaining random variables on the right hand sides of (32) - (34) are 
defined on a different probabihty space ^2,^2) , so that the random 
fields defined by the series are defined on the product probability space 
{n,T,¥) = (Oi X 02,-Fi X J"2,Pi X P2). Thus, for every fixed U2 G ^2, 
the conditional random field X((u;i, ^2)), t^^i G ^^i, is a centered Gaussian 
random field. We now apply to this random field Corollary 11.3.2 in [2]. 

Firstly we check the condition on the incremental variance of the second 
order partial derivatives there. In obvious notation, for every i,j = 1, . . . , d, 
and t, s e M, 

El [{X,j{t) - Xij{s)f} 
00 

= E i9^jit + Hk) - g^jis + Hu)f 1 {iP{VkMHk)Tk < 1) . 
fc=i 

Bounding the Holder constant of a function by the largest value of its partial 
derivatives, as in the proof of Lemma 6.4, we obtain 

El [{X,,{t) - X,,{s)f} 

00 

< (f\\t - sf J2 VkT^giHk)! {'^{VMHk)Tk < 1) , 

k=l 

and, hence, the incremental variance condition will follow once we check that 
the infinite sum above converges. For this, we need to check (see [14]) that 



^°°E2{min [1, V^Tg{Hfl{'iP{V)ip{H)x<l)'^} dx 



< 00. 



(The random variables without a subscript represent generic members of the 
appropriate sequences.) By the definition of the derivatives ip and ip, this 
reduces to checking that 



/ / min 1, y'^Tg{s)'^ dsdy < 00 



which is an elementary consequence of the integrability assumptions imposed 
on Tg in Condition 6.3, and of the assumptions imposed on the Levy measure 
Pq in Condition 6.2. 

It remains to check that the joint distribution under Pi of the random 
vectors of partial derivatives (Xj, Xij) is non-degenerate for P2-almost ev- 
ery a;2- This, however, follows from the representation (33) and (34) of the 
derivatives and the part of Condition 6.3 that rules out the possibility that 
the derivatives of the kernel g belong to a lower dimensional subspace. □ 
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6.2. Examples: How the shape of the kernel can affect the geometry of 
excursion sets. In the examples below Condition 6.2 is a standing assump- 
tion, and will not be mentioned explicitly. Our first example is of an infinitely 
divisible moving average random field whose high level excursion sets have 
a similar geometric structure to those of Gaussian random fields. 

Example 6.6. Let g be a non-negative kernel satisfying Condition 6.3, 
that is also rotationally invariant and radially decreasing. Formally, g{t) = 
gr{\\t\\) for some non-negative and decreasing function g,. on [0,oo). An ex- 
ample is the Gaussian kernel g{t) = exp{— a||t|p}, a > 0. It is trivial to check 
that the restrictions in Condition 6.3 on the various partial derivatives of g 
hold in this case. 

Corollary 5.1 tells us that the Euler characteristic of the excursion set 
over a high level, given that the level is exceeded, is asymptotically that of 
the field 

(35) {sup g{s + r))~^g{s + t), tGid, 

with a randomly chosen s € M"' and over a random level Va- 

The assumption of rotational invariance and radial monotonicity on the 
kernel g implies that, in this case, the excursion set of the random field is 
the intersection of a Euclidian ball centered at the point —s and the cube 
Ifi- This is a convex set and, hence, has Euler characteristic equal to 1, 
regardless of the point s € M*^ or the random level Va- 
in this case the limiting conditional distribution of the Euler characteristic 
is degenerate at the point 1. Furthermore, the excursion set has, with high 
probability, a "ball-like shape", as is the case for smooth Gaussian random 
fields. 

In spite of the 'Gaussian-like' conclusion in the previous example, it is easy 
to modify it to make the high level excursion sets of an infinitely divisible 
random field behave quite differently. Here is a simple example. 

Example 6.7. We modify the kernel g of the previous example by 
adding to it oscillations, while preserving its smoothness and integrability 
properties. For example, take 

g{t) = (1 + cos(0,t))e-''ll*ll', t € 

for fixed 6* G M"'. 

Then, depending on the random choice of the point s in (35), the structure 
of the excursion sets in could be quite varied, as it depends on the shape of 
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g in the translated cube . Thus, depending on the random level Va, the 
shape of the excursion set may be quite different from a ball-like shape. In 
particular, its Euler characteristic will have a non-degenerate distribution. 

6.3. The bottom line. The bottom line, of course, is that the shape of 
the excursion sets is determined, to a large extent, by the shape of the 
kernel in the integral representation of the random field or, alternatively, 
by the geometric properties of the functions on which the Levy measure of 
the random field is supported. By choosing appropriate parameters for the 
random field one can generate quite different distributions for the Euler and 
other geometric characteristics of high level excursion sets. 

Our hope is that this fact will generate greater flexibility in applications, 
allowing the practitioner to choose models with pre-determined excursion 
set shapes. 

Furthermore, the description of the limiting conditional distribution (and 
not only the expected value) of the numbers of critical points and so the 
Euler characteristic should allow one to devise better statistical tests based 
on the observed excursion sets. Precisely how this should be done, however, 
is something we shall leave for the future. 

7. Proof of Theorem 4.1 for the case M = [0, l]*^. The proof is 
rather long and rather technical, although the basic idea is not difficult. 

The basic idea, which is common to many proofs involving infinitely di- 
visible random fields X, is to write X as a sum of two parts, one which 
tends to be large and one which is made up of smaller perturbations. The 
large part, which, distributionally, behaves as a Poisson sum of deterministic 
functions with random multipliers, is comparatively simple to handle, and 
it is this part that actually accounts for the limit in Theorem 4.1. One then 
needs to show that the small pertubations can be ignored in the u ^ oo 
limit. In the argument that follows this is somewhat more difficult than is 
usually the case, since even if the small part is small in magnitude it can, 
in principle, have a major effect on variables such as the number of critical 
points of the sum. (Think of a constant function, f{t) = A, to which is added 
g{t) = €cos{{9,t)). No matter how large A might be, nor how small e might 
be, the critical points of f + g are determined by g, not /.) 

Due to the length of the ensuing proof, we shall do our best to signpost 
it as it progresses. 

(i) Some notation for the parameter space and for critical points. As men- 
tioned earlier, in this section we shall take as our parameter space the cube 
Id. The first step is to develop notation for describing its stratification. 
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Let Jk be the cohection of the 2^~^['^) faces of Id of dimension k, k = 
0, . . . ,d, and let J = [Jk •Jk- For each face J ^ Jk there is a corresponding 
set cr(J) C of cardinahty k and a sequence e(J) G {— l,!}'^^'^) 

such that 

J = [t = {ti, . . . ,td) e h tj = ej if j ^ a{J) and < < 1 if j G o-(J)}. 

Let g be a function on an open set M containing Id- For J ^ Jk and 
i = 0, 1, . . . , A:, let Cg( J; i) be the set of points i G J satisfying the following 
two conditions. 

(36) ^(t) = for each j G o-(J) , 

Otj 

(37) the matrix | ^ | has non-zero determinant 

\dtmdtn 

and its index is equal to k — i. 

Now define Ng(J; i) and Ng{J; i : u) in terms of Cg(J; z) as in Section 4. 

('iij Splitting X into large and small components. 

By Assumption 3.1, X and its first and second order partial derivatives 
are a.s. bounded on M, and, by (11), the Levy measure of X is concentrated 
on functions. Defining 

Sl = [g e C"^ : ma^ sup \g{t)\, jup \gi{t)\, ^ sup \gij{t)\ > l}, 

t£M teM ,i=l,...,d teM,i,j=l,...,d 

the sample boundedness of X, along with (14) and general properties of 
Levy measures on Banach spaces (e.g. [7]) imply that 



< oo. 



(38) d = Xx{Sl} 

We are now ready to decompose the infinitely divisible random field X 
into a sum of two independent infinitely divisible components by writing 

(39) X{t)=X^{t) + Y{t), teM, 

where X^ is a compound Poisson random field with characteristic functions, 
which, for /c > 1, ii, . . . , tfc G M, and real numbers 71, ... , 7fc, are given by 

(40) 

E I exp|i^7jXL(tj)| I = expj^ (^exp{i ^ 7jx(tj)} - l) Xxidx)^. 
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The second, or 'residual', component Y has characteristic functions 
E |exp{i^7jy(tj)}| = exp|-Q(7i,... ,7fc) 

+ / ~ {ew{i^ljx{tj)} - 1 - Ax((ix) +iLi(7i, 

JM. \Sl j^i j^i 

where we are using the notation of (13), and 

^i(7i,--->7fc) =^(7i,--->7fc) - / 

We shall ultimately show that the limiting behaviour of the critical points 
of X depends only on the component X^, so we study it first. 

(Hi) A limit theorem for the criticial points of X^ . We start by noting that 
it follows from the form of the characteristic function (40) and the definition 
(11) that X^ can, in law, be written as 

AT 

(41) X^(t)= 5^X^/(5fc,i), 

m=l 

where is a Poisson random variable with mean 6 given by (38), indepen- 
dent of an i.i.d. sequence of random pairs {{Xm, Sm)), m = 1, 2, . . .) taking 
values in (M \ {0}) x 5 with the common law 6~^F restricted to the set 

{(s,x) e (M\ {0}) X S : sup|x/(s;t)| > l} . 

tGA/ 

Recall that F is the Levy measure of the infinitely divisible random measure 
M in (4). 

Since the sum in (41) is a.s. finite, and the kernel / has bounded 
sections /(s; •) for all s € 5, it follows that X^ is bounded and on M. 

We now decompose the compound Poisson term X^ itself into a sum of 
two independent pieces, the stochastically larger of which will be responsible 
for the limiting behaviour of the critical points of X. For u > and 1/2 < 
P < 1 define the sequence of independent events 

Am{u) = ^max sup \Xmf{Sm]t)\, ^SUp \XrrJi{Sm]t)\, 



.7/ 



teM teM,i=l,...,d 



sup \Xjnfij{Sni;t)\ 
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and write 

N N 

(42) x^it) = + E^-/('5-;*n^™H= 

m=l m=l 

In Lemma 7.1 we shall show that X^^'"^^ and its partial derivatives have 
suprema the tail probabilities of which decay faster than the function H, 
and so are unlikely to affect the critical points of X. We shall return to this 
point later. 

Now, however, we shall concentrate on the critical points over high levels 
of X(^'i). Define two new events 

N N 

(43) Bi{u) = { Y l(^mN) = l}, B2{u) = [Y HAM) > 2} , 

m=l m=l 

The first of these occurs when there is a single large term in the Poisson 
sum (41), the second when there are more. On the event Bi{u) we define 
the random variable K{u) to be the index of large term, and otherwise allow 
it to be arbitrarily 

In the notation of Section 4 in general and Theorem 4.1 in particular, it 
follows that, on the event Bi{u), the following representation holds for the 
numbers of the critical points of X^^'^^ over the level u. For k = 0,1, . . . ,d, 
a face J ^ Jk and i = 0, 1, . . . , fc, 

Nj^(L,i){J;i : u) 

Ci{J\SK{u)) 

= H^K{u) > 0) Y ^{^K{u)f {SK{uy,ti{J;i; Sk{u))) > uj 

1=0 

Ck-i{J;SK{u)) 

+1(^XH<0) Y '^{XK(u)f{SK(uy,'ti{J;k-i;SK{u))) > ij) ■ 

1=0 

Therefore, for any number r = 1, 2, . . ., on the event Bi{u), we have Nj^(l,i) {J; i ■ 
u) > r if, and only if, 

Xk{u) > {f[r]'''^\SK(u)))~^U or Xk(^^) < -{fl;!^'''~'''~\SK(u)))~^U. 

We conclude that for any numbers n{J;i) = 1,2,..., for all J € J'k, and for 
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all /c = 0, 1, . . . , d and i = 0, 1, . . . , A;, 

(44) p{{7V^(L,i)(J;« : ii) > n(J;i) for allJ and i } n (n) } 

Write £m for the union of sets (— oo, — max] and [max, oo), where the 'max' 
come from the preceding lines with K{u) replaced by m. Then 



(45) p{{Xx(,) Gf^(„)}ni?i(u)} 



N 

U (Am{u) n fl Am, {uf n {Xm € £m}) 

m=l miy^m 



n=0 ' mi=2 

= er{Ai{u) n{Xie£i}}- p{{x^(„) g £Kiu)} n i?2( 

Applying this to the right hand side of (44) and using part (iii) of Lemma 
7.1 yields 

(46) p{{iVx(i„i)(^;« : u) > n{J;i) for ah J and i] nSi(n)} 

— <5small(^) , 



where Q sma\\{u) / H{u) ^ as n ^ oo. 

Assume for the moment that all the n( J; i) are strictly positive. Since the 
parameter (5 in the definition of the event Ai{u) is less than 1, it follows 
that, as u — > cx), 

p|^i(n)n(Xi> max {fl:'^n^l{Si)y^u\\ 

~p(Xi> max (/r'!'t^^l('Si))~^u| 



6* Js\' '^J&Jk,k-- 



max (/rl,/r\>i (-s)) ii, oo) mids). 
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In the last step we used the law of Xi introduced after (41) and the de- 
composition (5) of the measure F, and in the middle one the asymptotic 
equivalence means that the two ratio of the two probabilities tends to 1 as 
u — * oo. Since a similar asymptotic expression can be written for the second 
term in the right hand side of (46), we obtain 



lim 

U—>QO 



''{{^xiL.i){J;i : u) > n{J;i) for ah J and i} H Bi{u)^ 



H{u) 



lim H(u)-^ 
+p[s; (-00, 



P\ s; 



max 

JeJk,k=0,l,...,d, i=0 



,,„,(4iS,!(»r'". 



oo 



max 

JeJk,k=0,l,...,d,i=0 



m{ds) 



provided the last limit exists. Applying (17) and Potter's bounds, as in 
Lemma 3.6, to justify an interchange of limit and integration, and noting 
Assumption 3.5 relating p, and we have 



lim 



^{^[^x'^L,i){J]i : u) > n{J;i) for all J and i} H Bi{u)^ 



w+(s) min ( fy/^j^viis))"' 

+'W-{s) ^ ^ ^ . „, . (/| 



{J;k-i:-),-.^a 
J6j"fc,fc=0,l,...,d,i=0,l,...,fc'" ^ 



m{ds) 



Ic 



Finally, since by part (iii) of Lemma 7.1, the event B2{u) has a probability 
of a smaller order, we conclude that 



(47) 



lim 

M— >00 



^^N-^(L,i) (J; i : u) > n{J; i) for all J and zj 



as well. 

In view of (18), we can rewrite this as 



(48) lim P<^ 1) [J] i : u) > n(J; i) for all J and i I sup Xt>u\ 



This will complete the proof of the theorem, at least for the case of strictly 
positive n( J; i), once we show that the lighter-tailed random fields Y of (39) 
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and X^^'^^ of (42) do not change the asymptotic distribution of the numbers 
of critical points of X. This will take us a while to show, and makes up the 
remainder of the proof. 

Before we do this, note that handling situations in which some or all of 
the numbers n{J; i) are zero is actually only an issue of semantics, once we 
recall our convention regarding the 0-th order statistic introduced prior to 
the statement of the theorem. For example, in the case when all the n{J; i) 
are zero, the event on the left hand side of (48) should be interpreted as 
stating that X^^'^^ has crossed the level u, given that it has done so. Not 
surprisingly, the resulting limit, and the right hand side, turn out to be 1. 
Similar reductions work when only some of the n( J; i) are zero. 

(iv) An outline of what remains to do. It follows from what we have done so 
far that 

(49) X(t) =X(^'i)(t) + X(^'2)(t) + y(t), t£M, 
or, equivalently, that 

(50) X^^'^\t)=X{t)-X^^'^\t)-Y{t), teM. 

What we plan to show is that when when either X or X^^'^^ reaches a high 
level u, then the lighter-tailed random fields Y and X^^'"^^ can be thought 
of a small pertubations, both in terms of their absolute values, and those 
of their first and second order partial derivatives . This will imply that the 
asymptotic conditional joint distributions of the number of the critical points 
of the random fields X and X^^'^^ are not affected by the lighter tailed fields 
and, hence, coincide. 

In fact, what we establish is that near every critical point of one of the 
random fields X and X^^'^^ there is a critical point, of the same index, of 
the other. Equation (49) allows us to do this in one direction, and (50) will 
give us the other direction. The two equations are of the same type, and 
the fact that the terms in the right hand side of (49) are independent, while 
the terms in the right hand side of (50) are not, will play no role in the 
argument. Therefore, we shall treat in detail only one of the two directions, 
and describe only briefly the additional steps needed for the other. The first 
steps in this program involve collecting some probabilistic bounds on the 
closeness of critical points and the behaviour of Hessians there. 

(v) Bounds on critical points and Hessians. We start by introducing a func- 
tion D : S ^ {0, oo] that describes what we think of as the degree of non- 
degeneracy of the critical points of an s-section of the kernel /. This includes 
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the minimal Euclidian distance between two distinct critical points of an s- 
section of the kernel / and the smallest absolute value of an eigenvalue of 
the Hessian matrices of the section evaluated at critical points. Specifically, 
starting with critical points, and recalling the definition of the ti{J;i;s) as 
the critical points of index i on the face J for the s-section of X, define 

Di{s) = min|||ti^(Ji;ii; s) - t/al-^a; «2; : Jj G Jk^, 

0<ki,k2<d, < ij < kj, < Ij < Ci^{Jj;s), j = 1, 2}, 

where the minimum is taken over distinct points. Furthermore, define 

D2{s) = min||A| : A is an eigenvalue of {fmn{s;ti{J;i; s))) ^^^^^^j^ ; 

J ^Jk, 0<k<d, 0<i<k, l<l <Ci{J;s)Y 

As usual, both minima are defined to be equal to +00 if taken over an empty 
set. 

Now set 

(51) Dis) = inm{Di{s), D2{s)). 

Note that, by Assumption 3.4, D is a strictly positive function, so that for 
any any ^-valued random variable W one has liinT—>o^{D{W) < r} = 0. 
Choose W to have the law Nyy given by 

(52) ^{s) = c,K(s) + w.{s)) sup \f{s,tr, seS, 

where c* is a normalising constant. That this is possible is a consequence of 
(17). For e > 0, choose tq > so small that ¥{D{W) < tq} < e. With the 
random variable K{u) as before. Lemma 7.2 gives us that 
(53) 

r{{D{SKiu)) < TO, sup,g,^ |X(^'i)(t)| > u}nBi{u)} 
hmsup— 777-^ < e, 

where Bi{u) was defined at (43) and indicates that there was only one 'large' 
component in the decomposition of X. 
Note that, since the event 



< sup |x(^'^)(t)| > u\ n Bi{u) 



'tela 
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is a subset of Bi{u), on this event X^^'^\t) = Xx[u)f {S^(.^y,t) for all t € 

M. Thus, again on this this event, since the supremum of this field over 
Id exceeds n, while the kernel / is uniformly bounded, we conclude that 



l-'^xCit)! > ^^/||/||oo- Therefore, on the event 

{D(5^(„)) > TO, sup > n} n Bi{u) 

teid 

the smallest eigenvalue length 

(54) -Dmin = min||A| : A is an eigenvalue of fx^^^^(t)) ; 

L V / m,n(ia{J) 



J G iTfc, < k < d, t is a critical point on j| 



satisfies L>min > (to/H/Hoo)-"- 

We now combine (53) with (47) as follows. Introduce the event Qr{u) 
that occurs whenever the minimal Euclidian distance between two distinct 
critical points of the random field (^X^^'^\t), t G 1^) is at least r > 0, while 
the smallest eigenvalue length of the Hessian evaluated at the critical points 
satisfies -Dmin > (■'"/||/||oo)^*- Thus we have 

p|{iV^(L,i)(J;^ : u) > n{J;i) VJ, fl 17^„(ii)| 

(55) liminf^ — ^ > h-cZ^e, 

^ ' u^oo H{u) - * ' 

where Ic is as in (47). We can, furthermore, 'sacrifice' another e in the right 
hand side of (55) to add to the event a requirement that the largest 

eigenvalue of the Hessian evaluated at the critical points, which we denote 
by -Djnax) satisfies Dmax < Mu for some positive M = M{e). This is possible 
because -Dmax is bounded from above by the largest absolute value of the 
elements of the Hessian, which we bound from above by Mu with a large 
enough M. For the same reason, we can also bound from above the largest 
value of \\VX^^'^\t)\\ over Id by Mu. 

Denoting the resulting event by Qriu), we obtain 

(56) lim mf > - 2c, s. 

Now note that since, as stated above, X^ is bounded and on M, and 
the same is true for X by Assumption 3.1, it follows that the 'remainder' 
Y in (39) is also a.s. bounded and C^. Furthermore, by construction, Y 
and its first and second order partial derivatives have Levy measures that 
are supported on uniformly bounded functions. Consequently, the tail of 
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their absolute suprema decays exponentially fast; see [5]. In particular, for 
i,j = 1, . . . ,d, 

u^oo H{u) u->co H{u) 

~\Y,,{t)\>u] 
= hm -— 

M->00 H [U) 

= 0. 

It follows from this, part (ii) of Lemma 7.1, and the regular variation of H, 
that there is a function l{u) | oo such that l[u)/u ^ as u — cxo and 



(57) lim 



[sup^^~|y(i)|>/(^.)} 



M->00 Hiu) 



= lim 

ti— >oo 




\Y,,{t)\>l{u)} 




H{u) 


= lim 




\X^^'{t)\>l{u)} 




H{u) 


= lim 

M— ►oo 




\xl:f{t)\>Ku)] 




H{u) 


= 0, 







for i,j = l,...,d, 

We now combine (56) and (57) in the following way. Let n^r\u) be the 
intersection of the event ^}r{u) with the complements of all 4 events whose 
probabilities are displayed in (57) and set 

ncr{u) = {Nx^L.i){J;i:{l + T2)u)>n{J;i)yj,i} n n^^{u). 

Then, given < ei < 1, and using the regular variation of H, we can find 
Ti , T2 > such that 

,58) l,„„nfEM > (l-.O/e. 

u^oo H{u) 

(vi) The (almost) end of the proof. Continuing with the above notation, we 
now claim that, on the event Qcriu), for u large enough so that 

(59) — > max , — , 

l{u) \ n T2J 

we also have 

(60) Nx{J;i:u)>n{J;i), J £ Jk, k = 0,1, . . . ,d, i = 0,1, . . . ,k . 
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Note that, once this is estabhshed, we shaU have 



hm inf 



F{Nx{J;i : u) > n{J]i), J <£j, < i < dim J} 



> {l-ei)h 



C) 



and, since this holds for ah < ei < 1, we also have 



(61) 



lim inf 




> Ic. 



H{u) 



Combining this with (18) gives Theorem 4.1, albeit with an inequality 
rather than an equality in (24). 

To obtain the opposite inequality assume that, to the contrary, for some 
numbers n( J; i). 



along some sequence j oo. 

Now proceed by repeating the steps performed above and, this time using 
(50) rather than (49), and so demonstrate the existence of a critical point of 
near each one of X. Thus (62) also holds with X replaced by X^^'^^), 



Since this contradicts (47), (62) cannot be true, we have the required lower 
bound, and the proof of Theorem 4.1 is complete, modulo the need to es- 
tablish the claim (60). 

(vii) Establishing (60) to finish the proof. In order to establish (60), we 
shall show that, on the event Q.cr{u)^ to every critical point above the level 
(1 + T2)u of the random field X^^'^^ we can associate a critical point above 
the level u oi X which is in the same face and of the same type. 

To this end, let to be a critical point above the level (1 + T2)u of X^^'^'^ 
that belongs to a face J G cTfc for some < k < d, and which is of the type i 
for some < i < k. Let (ei, . . . , e^) be an orthonormal basis of M*^ consisting 
of normalised eigenvectors of the Hessian matrix 



(62) 




VIZ. 



n 




(63) 
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and let Ai, . . . , be the corresponding eigenvalues. Note that, by the def- 
inition of the event ^lcr{u), we have |A„,| > (Ti/||/||oo)ii for n = 1, . . . ,/c. 
We naturally embed the vectors (ei, . . . , e/c) into the face J and make them 
d-dimensional vectors by appending to them the d — k fixed coordinates of 
the face J. (We shall continue to denote these vectors by (ei, . . . , efc).) Note 
that for small real numbers ei , . . . , we have 

k k 

(64) VX(^'i)(to + e.ej) = ^ e,A,e, + o(max(|ei|, . . . , \ek))) • 

i=i i=i 

In particular, the directional derivatives 

9f'^)(t) ^(Vx(^'i)(t),e,), j = l,...,fc, 

satisfy 

k 

(65) 9j^'^\to + J2^j^j) = ejAj +o(max(|ei|,...,|efc|)) . 

In what follows we shall work with a small positive number e > 0, placing 
more and more conditions on it as we progress, to clarify precisely how small 
it will need to be. As a first step, take e < ti/2, where ri is as in (58). 

Consider a A;-dimensional cube (which is a subset of the face J) defined 

by 

k 

Ce = {to + Y9jej,\ej\<e,j = l,...,k^, 

3=1 

along with its {k — l)-dimensional faces 

k 

Fn = {^0 + E = ±e, \0j\ < e, 1 < J < k, j ^ n] , 

i=i 

where n = 1, . . . , A;. It follows from (65) that, for e > small enough, n > 1, 
and, as above, M large enough, we have 

(66) 2Meu > 2e\\n\ > l^f '^)(t)| > ^ > 

for all t € Fl^^ n = 1, . . . ,k. The assumption that e be small enough now 
entails that (66) holds for all critical points and for all relevant n. Since 
the number of critical points is finite, this requirement is easy to satisfy. 
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Similarly, the continuity of the eigenvalues of a quadratic matrix in its com- 
ponents (see e.g. Section 7.2. and Corollary 2 in Section 7.4 of [6]) shows 
that, for all e > small enough, the eigenvalues of the matrix of the second 
order partial derivatives {Xmn{t)m,n£a{j) have all absolute values satisfying 
|A„| > (r/2||/|| oo)u for n — 1, . . . , fc and t E C^. Finally, we require that e 
be small enough that this lower bounds holds for all critical points to con- 
sidered above. In particular, this implies that the signs of these eigenvalues 
throughout are the same as those at the point to- 

Next, for a non-empty / C {1, . . . , fc} and p £ { — 1, 1}'"' consider the vector 

(67) x{I,p) = ^piei. 

Consider a point t that belongs to the (relative to the face J) boundary of 
the cube and, more specifically, belongs to the face of that cube defined 
by 

(68) ( n n ^' 

ieI,Pi=l i£l,Pi=-~l 

and to no other (k — l)-dimensional face of C^. Define a function h^^'^^ 

i=l 

This is a C^-function, and its gradient (within the face J) is given by 

V/.(^'^)(t) = 2X:<7f ''^(i)V#'^)(t) = 2;^#'^)(i)7^(^'i)(i)ef . 
j=i j=i 

Note also that for all / and p as above. 



if f ^ /. 



V5f'')(t0),rE(/,p)) 



In particular, we can write for any t belonging to the face of defined by 
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(68) , 

(69) (Vh^^'^\t), x{I,p)) 

= 2^5f''^(t)(V5p)(t),x(/,p)) 

i=l 

+2j:gi'^'\t) ((V#'^)(t) - V#'^)(to)) , xil,p)) . 

i=l 

It follows from (65) and (66) that 

(70) gi^'^\t) > for i G F+ if A„ > and for t G F~ if A„ < 0, 
gi^'^^^t) < for i G F+ if A„ < and for t G F~ if A„ > 0. 

Consquently, we can conclude, by (70) and (66), that the first term in the 
right hand side of (69) is negative and, more specifically, does not exceed 

-2Card(/)Z)^i„^^n < -{Ti/\\f\Uhu\ 

We can bound the absolute value of the second term in the right hand side 
of (69) from above by 

2fcE bf ''^(t)l • ||V#'^)(t) - Vgl'''\to)\\ < 2k'MV, 
1=1 

by the definition of the event i^criu)- This, obviously, indicates that, for 
e > small enough, 

(v/i(^'i)(0, x{I,P)) < -Ceu\ 

where C is a finite positive constant determined by the parameters in the 
event Qcriu). If 

gj{t) = {VX{t),ej), j = l,...,k, 

and we define 

(71) h{t)=j2{m{t))\ 

i=l 
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then, on the event Qcr{u), 

{Vh{t), x{I,p)) < -Ceu^ + kl{uf. 

Taking into account that l{u)/u — > as n ^ oo, where / is given by (57), 
we see that for u large enough it is possible to choose e > small enough 
such that 

(72) (V/i(t), x{I,p)) < 0, 

for any t belonging to the face of defined by (68). The final requirement 
on e is that (72) holds. 

Similarly, since by the definition of the event ^}cr{u), the first order partial 
derivatives of X^^'^^ and Y are bounded by l{u) = o(u) in absolute value, 
we have that (70) and (59) give us 

(73) gn{t) > for t G F+ if A„ > and for t G F' if A„ < , 
gn{t) < for t G F+ if A.„ < and for t G if A„ > 

as well. 

In order to complete the proof and establish (60) it is enough to prove 
that, on Qcriu), X has a critical point in the cube Cg. In fact, if such a 
critical point exists. Lemma 7.3 below implies that it will be above the level 
u and of the same type as to- Furthermore, these critical points of X will all 
be distinct. 

To establish the existence of this critical point, note that, by the continuity 
of \/X and the compactness of C^, there is a point ti in at which the 
norm of the vector function g{t) = (51 (t), . . . ,gkit)) achieves its minumum 
over Cg. We shall prove that, in fact, g{ti) = 0. By the linear independence 
of the basis vectors ei, . . . , e^, this will imply that gj{ti) = for j = 1, . . . , fc, 
and so ti is, indeed, a critical point. 

Suppose that, to the contrary, g{ti) ^ 0, and consider firstly the possi- 
bility that the point ti belongs to the (relative to the face J) interior of C^- 
Note that the Jacobian of the transformation g : — > M'^ is given by 

jg{t) = £n{t), 

where fi^t) = (^mn(i))m nGo-(j) Hessian of X, and £ is a k x k matrix 

with rows ei, . . . , e^. We have already established above that, on the event 
^Icriu), TC is non-degenerate throughut C^. Since the vectors ei, . . . , are 
linearly independent, we conclude that the matrix £ is non-degenerate as 
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well. Since the vector g{ti) does not vanish, it has a non- vanishing compo- 
nent. Without loss of generality, we can assume that 51 (ti) 7^ 0. Choose a 
vector 2; G M'^ for which 

Jg{h)x' = (1,0,..., 0)'. 

Then for 5 G M, with \5\ small, 
g{ti + 5x) = 



and so 

1 1 2 

||5f(ti + = 

< 

for 6 with \5\ small enough and such that dgi{ti) < 0. This contradicts the 
assumed minimality of 11(7(^1)11 and so we must have g{ti) = 0, as required, 
for this case. 

It remains to consider the case g{ti) 7^ 0, but the point ti belongs to the 
boundary of the cube C^- Let g{ti) belong to the face of the cube defined 
by (68). With the function h defined in (71), we have, for 5 > small, 

h{ti+5x{I,p)) = h{ti) +6{Vh{t),x{I,p)) + o{6) . 

By (72), this last expression is smaller than /i(ti) if 5 > is small enough. 
However, by the definition of the vector x{I,p), the point ti+5x{I,p) belongs 
to Ce for 5 > small. Once again, this contradicts the assumed minimality 
of ||g(ti)||. 

Thus we have established (60) and, therefore, (61), and so the Theorem, 
modulo the need to prove the following three lemmas. □ 

Lemma 7.1. The following three results hold: 
(i) The random fields X^^'^^ and X^^''^^ on the right hand side of the de- 
composition (42) are independent. 



g{ti) + 5Jg{ti)x^ + oi\S\) 
( gi{ti) + 5 \ 

92{tl) 



+ o{\6\), 



\ 9k{ti) J 



J2g^{ti)^ + 25gi{t,) + om 

k 
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(a) The random field X^^''^'^ has sample functions and satisfies 
^.^ ¥{snp^^~\X(L'^){t)\>u} ^ P{sup^^~|Xp)(t)|>n} 

u^oo H{u) u^co H{u) 

Ji^oo H{u) ' 

i,j = l,...,d, where H is the regularly varying function of Assumption 3.5. 
(Hi) The number of terms in the sum defining X^^'^^ satisfies 

hm — = 0. 

Proof. The claim (i) follows from the fact that a Poisson random mea- 
sure, when restricted to disjoint measurable sets, forms independent Poisson 
random measures on these sets (see e.g. [12]). Since the sum defining the 
random field X^^'"^^ is a.s. finite, the fact that it has sample functions in 
follows from Assumption 3.1. Furthermore, for e > 0, choose > so large 
that P{A^ > n^} < e. The above discussion implies that the number K{u, e) 
of the terms in the sum defining X^^'^^ in (42) that satisfy 

u 

sup \XmfiSm;t)\ > — 

is Poisson with the mean less or equal to 

f|(s,x) G (M\ {0}) X S : sup \xf{s\t)\ > u/(2?i,)} 

= \x{9: sup|5(t)| >n/(2n,)} - CH{u) 

teM 

as u ^ oo, where we have used (11) and Lemma 3.6. Therefore, for large u 
P{sup \X^^'^\t)\ >u} < F{N > nj + ¥{K{u, e) > 2} < e + CH{uf, 

and so 



P{sup 

limsup ^ ^ . . < e. 

Letting e — > completes the proof of the first limit in part (ii) of the lemma, 
and the other limits are established in the same way. Part (iii) of the lemma 
can be proven similarly. □ 
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Lemma 7.2. The random variables Dm in (51) satisfy 

f{{^(5a'(„)) < TO, sup,^,^ > n} n B^iu)} 
limsup — — < e , 

where is as in (52). 

Proof. We use a decomposition as in (45) to obtain 
¥{{D{SKiu)) < TO, sup|x(^'^)(t)| > u}nB,{u)} 

n 

mi =2 ^^^d 

Since on the event Am{uy one has sup^^^^ \Xmf{Sm:,t)\ < u^, it follows 
that the latter probability can be asymptotically bounded by 



e¥{D{Si) < TO, sup|Xi/(Si;t)| > u} 

tela 

= / / t{D{s) < To, \x\snp\f{s;t)\ > u) F{ds,dx) 
Js Jm.\{o} tela 

= I l{D{s)<To) \ [ l(\x\> u(svLp\f(s;t)\)'^) p(s;dx)] m(ds 
Js \M{0} ^ tela ' ) 



^H{u) / l{D{s)<To){w+{s) + w.{s))sup\fis;tTm{ds), 
Js tela 

where we used Assumptions 3.5, 3.4 and (17). The lemma now follows from 
the choice of tq. □ 

Lemma 7.3. Suppose that for every critical point to of the random field 
(X^'^{t), t S Id), the random field {X{t), t € Id) has, on the event Qcr{u), 
a critical point in the cube ■ Then the critical points of X in Id correspond 
to distinct critical points of X^'^ , are themselves distinct, are all above the 
level u, and each of them is of the same type as the corresponding critical 
point of X^'^ . 

Proof. The fact that the critical points of X corresponding to distinct 
critical points of the field X^'^ are all distinct follows from the lower bound 
on the distance between two distinct critical points of X^'^ in the definition 
of the event i^criu) and the choice of e. The fact that all the critical points 
are above the level u follows from the lower bounds on the values of X^'^ at 
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its critical points in the definition of 0,criu) and, once again, the choice of 
e. It remains, therefore, to prove that a critical point in the cube of X is 
of the same type as the critical point to of X^'^. 

To this end, note that the absolute values of the eigenvalues of the matrix 
of the second order partial derivatives {X^^{t) + ^n(*))„^ are, on 

the event ilcriu), bounded from above by 2kl{u). Using continuity of the 
eigenvalues of a quadratic matrix in its components (see, once again. Section 
7.2 and Corollary 2 in Section 7.4 of [6]) we see that the Euclidian distance 
between an eigenvalue of (^Xmn{t)) ^ nGo-(j) ^^'^ corresponding eigenvalue 
of (-'^mn(0)m necr(j) ^® bounded from above by 2k l{u). Using the choice of e 
then shows that the numbers of the negative eigenvalues of the two Hessians 
are identical, as required. □ 

8. Proof of Theorem 4.1 for the general case. If you carefully 
followed the proof of Theorem 4.1 for the case M = 1^, then you will have 
noticed that there were two main components to the arguments. In the first, 
global bounds on suprema of random fields, derivatives, etc, played a crucial 
role. These arguments are no different for cubes than they are for other 
compact sets, including compact stratified manifolds, and so throughout we 
could have worked with general M rather than the special case M = Id 
when handling these arguments. 

The second component relied on showing that the critical points above 
high levels of the random fields X and X^'^ were in one-one correspondence, 
in terms of their (approximate) positions, their (approximate) heights, and 
their type. The arguments here were of a purely local nature, and took 
into account, for example, on what face of the critical points occured. 
Transfering these arguments to the case of M a stratified manifold is not 
trivial, but it is also not too hard. The main step lies in obtaining an analogue 
to the linearisation (64) of VX^'^ (the random field itself is defined precisely 
as in the simpler case), as well as an analogue of the Hessian (63). 

The natural place to do this, of course, is on the tangent spaces at the crit- 
ical points of X^'^. In small neighbourhoods of each such critical point, both 
X and X^'^ can be pushed forward to the tangent space via the exponential 
map. Once this is done, we are are back in a simple Euclidean setting and 
can argue as before, but now with respect to the push forwards. Pulling the 
results back to M itself, again via the exponential map, is straightfoward, 
as long as the mapping is smooth enough. However, looking again at the 
case M = Id, it is easy to see that all that 'smooth enough' really requires is 
a universal bound on the second order partial derivatives of X^'^ and that 
the pullbacks can be done in a uniform fashion. It is to ensure this this is 
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indeed the case that the assumption that the component manifolds djM of 
M have bounded curvature, made aheady in Section 2, was introduced. 

We leave the details to the interested reader. If you are uninterested, or 
unprepared to take our word for the fact that all that remains is, in essence, 
to repeat the above argument with heavier notation, you have the choice of 
believing only the case M = Id- There is enough new in this case to make 
the basic results of this paper interesting and useful. 
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